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Analyzing the spectrum of general, non-hermitian Dirac operators 
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We discuss the computational problems when analyzing general, non-hermitian matrices and in particular the 
un-modified Wilson lattice Dirac operator. We report on our experiences with the Implicitly Restarted Arnoldi 
Method. The eigenstates of the Wilson-Dirac operator which have real eigenvalues and correspond to zero modes 
in the continuum are analyzed by correlating the size of the eigenvalues with the chirality of the eigenstates. 
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1. THE PROBLEM 

Analyzing the eigensystem of large matrices 
such as the Dirac operator in 4-D lattice field 
theory is a quite demanding numerical task. Fur- 
thermore the lattice Dirac operator in Wilson for- 
mulation is neither hermitian nor anti-hermitian 
which poses an additional difficulty compared to 
a case with symmetry. 

The Wilson-Dirac operator consists of a triv- 
ial part (proportional to the unit matrix) which 
comes from mass- and Wilson terms and a non- 
trivial part K. Here we analyze the clover- 
improved jt) Wilson-Dirac operator and K con- 
sists of two parts K = Q — c sw C where Q is 
the standard Wilson hopping matrix and C the 
clover term. The matrix K is neither hermitian 
nor anti-hermitian, but hermitian conjugation is 
implemented through the similarity transforma- 
tion 
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The relation (Q) allows to define the hermitian 
modification 75 -?T and several numerical studies 
of this modified matrix can be found in the lit- 
erature. For a deeper understanding of the lat- 
tice Dirac operator however one would also like 
to study the original problem. 
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The two main techniques to solve large scale 
eigenvalue problems are Krylov subspace meth- 
ods: The Lanczos algorithm, which transforms 
the original matrix to tridiagonal form, or, in 
the case of non-hermitian matrices, the Arnoldi 
method, which transforms to Hessenberg form. 
Since typically the Hessenberg form is dense it is 
necessary to compress the relevant information in 
a Krylov subspace of small dimension. This can 
be done by using the Implicitly Restarted Arnoldi 
Method (IRAM) §]. IRAM has several favorable 
features: The numerical precision is comparable 
to dense methods, IRAM is able to find degen- 
eracies and it gives both eigenvalues and eigen- 
vectors. The software package ARPACK with 
an efficient and reliable implementation of the 
IRAM was developed by Sorensen and collabo- 
rators ||. The method turns out to be particu- 
larly useful when analyzing the edges of the spec- 
trum. It can e.g. be set up such that the eigen- 
values with largest real parts are computed first. 
This gives the eigenvalues in the physical branch 
of the spectrum. In Fig. [t] we show the physi- 
cal branch of the spectrum in the complex plane 
for a thermalized SU(2) gauge field configuration 
(12 4 , (3 = 2.4, c sw = 1.4) and compare it to its 
cooled counterpart (c sw = 1). Both gauge field 
configurations were taken from B. 
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Figure 1. Physical branch of the spectrum of 
K for a thermalized configuration and its cooled 
counterpart. 

As is obvious from the plot cooling leads to a 
strong shift of the eigenvalues and orders them 
along a single curve (ellipse). What is also in- 
teresting to note is the fact that the thermalized 
configuration has considerably more real eigenval- 
ues in the physical branch than its cooled coun- 
terpart. In m it was shown, that these ad- 
ditional real eigenvalues are partly due to the 
clover term. This proliferation of real eigenval- 
ues makes a probabilistic interpretation of the 
Atiyah-Singer index theorem for thermalized con- 
figurations problematic, as will be discussed in 
the next section. 

For a complete understanding of the properties 
of the lattice Dirac operator it is also interesting 
to analyze eigenvalues in the interior of the spec- 
trum, in particular on the real axis. However, 
as with all Krylov subspace methods, the conver- 
gence of the interior eigenvalues is very slow and 
one has to use some kind of spectral precondi- 
tioning. Traditionally, to find the eigenvalues in 
the vicinity of some point z in the complex plane 
a shift of the origin is combined with inversion, 
i.e. the eigenvalues of (K — are searched 

for. This requires solving a large linear system at 
every iteration which makes this approach quite 
inefficient. However, for matrices which satisfy 
(|l|) a kind of sophisticated preconditioning can 
be used ||. The shift of the origin is combined 



with a multiplication with 75 to transform the 
original non-hermitian to a hcrmitian eigenvalue 
problem. Although the explicit relation between 
eigenvalues of original and transformed matrix is 
not known the zero eigenvalues and their eigen- 
vectors are the same. By shifting the origin along 
the real axis, it is possible to identify real eigen- 
values. In particular one traces the flow of low 
lying eigenvalues of H{p) = 75 [AT — pi] as a func- 
tion of p. Whenever the flow crosses zero at 
some p — r, then this r is a real eigenvalue of 
K. The chirality of the mode can be computed 
from the slope of the eigenvalue flow at the cross- 
ing. Even for this hermitian problem we found 
the IRAM method superior (numerical stability, 
precision and efficiency) to our implementation of 
the accelerated conjugate gradient algorithm Q 
which is usually used for this purpose. Further 
gain in efficiency can be achieved by using Cheby- 
shev polynomials to project out the desired part 
of the spectra. 

3. SOME RESULTS ON THE REAL 
PART OF THE SPECTRUM 

It was already remarked that the real eigenval- 
ues of the lattice Dirac operator are of particu- 
lar interest. This is a consequence of (Q) which 
implies that tp^j^ip is different from only for 
eigenvectors tp of the Dirac operator which have 
real eigenvalues. Thus a comparison of the chi- 
ral properties of the eigenstates implies that only 
eigenvectors with real eigenvalues can play the 
role of the zero modes in the continuum [^|J^]. 
Based on this interpretation of the eigenvectors 
with real eigenvalues as the lattice zero modes 
one can formulate a lattice version of the Atiyah- 
Singer index theorem: v\U] = R- —R+. Here R+ 
and R- are the numbers of real eigenvalues in the 
physical branch of the spectrum with positive and 
negative chirality. The chirality is defined as the 
sign of the pseudoscalar matrix element of the 
corresponding eigenvectors and u[U] denotes the 
topological charge of the gauge field. 

The lattice version of the index theorem is 
expected to hold for sufficiently smooth gauge 
field configurations. This was demonstrated for 
smooth SU(2) and SU(3) background configura- 
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tions U and for fully quantized QED 2 (l0|. For 
thcrmalized SU(2) and SU(3) configurations the 
situation is less clea r There is no one-to- 

one corespondence [^[between the value of the 
topological charge assigned using the index theo- 
rem and the value from improved cooling [Q] and 
only averaging over larger samples or studying 
density functions of eigenvalues allows a prob- 
abilistic interpretation of the index theorem. 

HereQ we correlate the position of the real 
eigenvalues with the size of the pseudoscalar ma- 
trix element tp^j^ijj for the corresponding eigen- 
vectors ip. For smooth gauge field configurations 
the real eigenvalues in the physical branch of K 
are in the vicinity of 4 on the real axis and the 
pseudoscalar matrix elements are in the vicin- 
ity of ±1, the sign given by the chirality of the 
mode. For thermalized configurations the abso- 
lute value of the pseudoscalar matrix elements be- 
comes smaller and additional pairs of real eigen- 
values with opposite chirality are created || . We 
check if the artificial real eigenvalues from the 
clover term have smaller pseudoscalar matrix el- 
ements and if this can be used to discriminate 
between the would-be zero modes on the lat- 
tice and the additional real modes. In Fig. || 
we show a scatter plot containing the eigenval- 
ues between 1.8 and 4 on the real axis using a 
sample of 10 thermalized gauge field configura- 
tions (12 4 ,/3 = 2.4, c sw = 1.4). The horizontal 
axis gives their position on the real axis and on 
the vertical axis we plot the pseudoscalar matrix 
elements of the corresponding eigenvectors. 

If there was a clear separation of the lattice 
would-be zero modes from artificial real modes 
(due to strongly fluctuating fields and the clover 
term) one would find two well isolated islands of 
points in the upper and lower right corners of our 
plot. However, as is clear from the figure no such 
separated islands exist and we have to close with 
the conclusion of |5| that for 12 4 ,/3 = 2.4, c sw = 
1.4 there is no proper separation between physical 
zero modes and artifacts. 




Figure 2. Correlation of the real eigenvalues with 
the corresponding pseudoscalar matrix elements. 
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